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[57] ABSTRACT 

Periodic gain adjustment in plants of irreducible order, n, or 
for equalization of communications channels is effected in 
such a way that the plant (system) appears to be minimum 
phase by choosing a horizon time N>n of liftings in periodic 
input and output windows p„ and p >; , respectively, where N 
is an integer chosen to define the extent (length) of each of 
the windows p w and p y , and n is the order of an irreducible 
input/output plant. The plant may be an electrical, mechani- 
cal or chemical system, in which case output tracking (OT) 
is carried out for feedback control or a communication 
channel, in which case input tracking (IT) is carried out. 
Conditions for OT are distinct from IT in terms of zero 
annihilation, namely H ? H + V =I for OT and H S + H=I of IT, 
where the OT conditions are intended for gain adjustments 
in the control system, and IT conditions are intended for 
equalization for communication channels. 
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EXTENDED HORIZON LIFTINGS FOR 
PERIODIC GAIN ADJUSTMENTS IN 
CONTROL SYSTEMS, AND FOR 
EQUALIZATION OF COMMUNICATION 

CHANNELS 5 

ORIGIN OF INVENTION 

The invention described herein was made in the perfor- 
mance of work under a NASA contract, and is subject to the 10 
provisions of Public Law 96-517 (35 USC 202) in which the 
contractor has elected not to retain title, 

TECHNICAL FIELD 15 

The invention relates to providing periodic gain adjust- 
ment in control design as a way of modulating control 
energy to higher frequencies and making the harmonics add 
in such a way that a nonminimum-phase system appears to 
be minimum phase using a technique based on the notion of 20 
a mathematical “lifting” in which a serial-to-parallel con- 
version is performed on the control input and output signals, 
and mappings are considered between the vectorized quan- 
tities. This lifting technique is also applicable in communi- 
cations for equalization of a nonminimum phase channel. 25 

BACKGROUND OF THE INVENTION 

In the prior art, nonminimum-phase systems could only 
be controlled by using low bandwidth, sluggish controllers. 30 
In layman’s terms, nonminimum-phase systems move in a 
direction initially opposite to the direction they are pushed. 
Such systems arise when, for example, one attempts to 
dampen vibrations on a flexible structure in optical instru- 
ments by acting at a location remote from the flexible 35 
structure in the case of a “noncolocated” control system. 

Nonminimum-phase systems arise in many aerospace 
applications of noncolocated control, such as spacecraft 
slewing, instrument pointing/tracking, flexible robotics, 4Q 
acoustic systems, chemical process control systems, com- 
pact disk controllers, floppy disk controllers and even real- 
time equalization of nonminimum-phase communication 
systems. An object of this invention is to achieve high 
performance control for the class of nonminimum-phase 45 
systems that encompasses these and other applications 
which are referred to genetically as “plants.” 

A Zero Annihilation Periodic (ZAP) control law was 
introduced by the present inventor, David S. Bayard, in 
publications titled “Globally stable adaptive periodic con- 50 
trol,” Jet Propulsion Laboratory, Internal Document JPL 
D-9448, February 3, 1992 and “Zero annihilation methods 
for direct adaptive control of nonminimum-phase systems,” 
Proc. Seventh Yale Workshop on Adaptive and Learning 
Systems, Yale University, May 1992, for controlling non- 55 
minimum-phase systems using stable inversion of the plant. 

The general ZAP approach is based on the notion of a 
mathematical “lifting” performed on the plant input and 
output signals with mapping between the vectorized input 
and output quantities. 60 

A key property of the ZAP liftings which makes them so 
useful is that the transmission zeros of the lifted and vec- 
torized plant input signal are annihilated (i.e., placed to the 
origin). This zero-annihilation (ZA) property allows the 
vectorized plant signals to be stably inverted using standard 65 
control methods. The result is important to many areas of 
control, communications, and signal processing where a 


2 

stable plant inverse is often desired but not possible due to 
nonminimum-phase restrictions. 

In a paper by R. Lozano-Leal, “Robust adaptive regula- 
tion without persistent excitation,” IEEE Trans. Automatic 
Control, Vol. 34, pp. 1260-1267, December 1989, a multi- 
rate sampling method is presented which allows stable 
inversion of any linear time-invariant finite- order plant 
signal. Not surprisingly Lozano’s lifting has also been 
applied to developing stable adaptive control algorithms for 
nonminimum phase systems. See also R. Lozano et al., 
“Singularly-free adaptive pole placement using periodic 
controllers,” IEEE Trans. Automatic Control, Vol. 38, pp. 
104-108, January 1993. 

Lifting in this invention is an extension of the prior art 
represented by Bayard, supra, and Lozano-Leal, supra, to 
“extended horizons,” i.e., to extended periodic windows 
over which liftings (samples) are taken in adaptive control- 
lers of the type described by Lozano et al., supra, for 
example. Extended horizon liftings are crucial for control 
gain reduction in order to allow practical implementations of 
control in nonminimum-phase systems. The distinction 
between extended horizon liftings of the present invention 
and the prior art will now be described. 

Since Lozano-Leal’ s liftings are distinct from the present 
invention in that it utilizes a horizon size of N=2n, where n 
is the plant order and N is the length of the window for the 
liftings, and the number of elements in the system input 
equals those in the system output (i.e., a u =a > ), it will be 
denoted as the “2n-lifting ,” while in the present invention 
N>n and preferably N^2n-1 in a large class of systems in 
which a„>a y for output tracking (mapping), such as in 
feedback control systems, and G u <c y for input tracking, such 
as in communications systems, and which enjoy the same 
zero annihilation properties as N=2n. Unlike the special case 
of Lozano-Leal’ s and other prior art liftings, the present 
invention allows the use of all extended horizon liftings with 
N greater than n and preferably equal to or greater than 2n-l 
with cs u >G y or output tracking (OT) and G u «J y for input 
tracking (IT). An important consequence is that such 
extended horizon liftings lead to plant-inverse controllers 
with significantly reduced control gains. This overcomes a 
problem associated with Lozano-Leal’ s lifting, where N=2n 
and which has prevented its use in many applications 
of practical interest that the present invention will reach. To 
illustrate the present invention, a simulation example is 
provided below in which the peak control requirement is 
reduced by four orders of magnitude using extended horizon 
liftings. 

It will also be shown that as a dual result, a related class 
of extended horizon liftings enables equalization of non- 
minimum-phase channels in communication systems. This 
overcomes the standard problem of inverting the channel in 
a stable fashion. In this invention, the extended horizon 
property allows channel inversion by least squares estima- 
tion, which provides smoothing in the case of noise. 

As background information, consider a plant input/output 
model, 

Kz^)yi = BU~ J )u t (l) 

A{z~ l ) = 1 + £ a*-*; B(z~ ] ) - £ bur* 

£=1 i'=l 

where the polynomials A and B are assumed to be relatively 
prime. It is assumed that b^O, so that the polynomial B can 
be factored uniquely into the form B(z'~ 1 )=z -rf b 1 B(z _1 ) 
where B(z -1 ) is a monic polynomial and d-1 is the plant 
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delay. The choice d=T is for convenience only and is not a 
fundamental restriction. In the case that d^l, all subsequent 
expressions can be appropriately modified without loss of 
generality. 

Choose some horizon time N>n. The system of Eq. (1) is 5 
iterated to give the following system of linear equations, 


Yi^AJikhAjYik-D+By U(k)+B 2 U(k -1) (2) 


where, 


10 


1 I 

vjw+1 


• « 

UkN 

ykN+2 


UkN+1 


; m = 


yicN+N 

k J 


MJ W+N-\ 

a a 


Y(k)=plant output 
u(k)=plant input 

A flower triangular Toeplitz, with first column [0,-ai, . . . 

... , 0 ] 

A 2 =upper triangular Toeplitz, with first row [0,-ap . . . 

-a„A ... ,0] 

B flower triangular Toeplitz, with first column [b 1 ,b 2 , . . . 
,b n A • ■ Of 

B 2 =upper triangular Toeplitz, with first row [O, . . . , O, b n , 
• ■ . , b 2 ] 

Example 1 Let n=3 and N=4. Then Eq. (2) becomes. 
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It is convenient to combine terms involving Y(k) in eq. (2) 
and rearrange to give the following lifting of P. Albertos, 
“Block multirate input-output model for sampled- data con- 55 
trol systems,” IEEE Trans. Automatic Control, Vol. 35, No. 

9, pp. 1085-1088, September 1990. 

Albertos’ Lifting: 


Y(k)^AY(k-\)+HU(k)+BU(k-l ) (4) 

where, 


A=(I-A t r l A 2 

(5a) 

H=a~A t r i B t 

(5b) 

B=(I-A l )-'B 1 

(5c) 


60 


65 


4 

It is noted that since A x is lower triangular with zeros on the 
diagonal, the quantity (I-A^ is always invertible. Hence the 
quantities in Eq. (5) always exist. 

Polynomial A is divided into B to give impulse response 
sequence {hj, 




OQ 


= Z hiZ-t 
r=l 


(6) 


The Markov parameter sequence {hj is not assumed to be 
convergent (i.e., the system may be unstable). Using the 
Toeplitz structure of and Bj and relation, (Eq. (6)), it can 
be shown that the matrix H in Eqs. (4) and (5) can be written 
in terms of the impulse response parameters, 


h\ o 

hi hi 


hff . . 



(7) 


This is the desired expression for H, i.e., H=lower triangular 
Toeplitz, with first column [hj,h 2 , . . . , h H ] T . Since the delay 
is unity by assumption (i.e., d=l), the matrix H has a nonzero 
diagonal (i.e., h^), and is always invertible. 


STATEMENT OF THE INVENTION 

In accordance with the present invention, an improved 
method for providing periodic gain adjustment in plants of 
irreducible order, such as for control of electrical, mechani- 
cal or chemical systems or for equalization of communica- 
tion channels, comprises the choosing of a horizon time N>n 
of liftings in periodic input and output signal windows p„ 
and p y , respectively, where N is an integer chosen to define 
the extent (length) of each of the windows p w and p v , and n 
is the order of the irreducible input/output plant model given 
by an equation of the form 


A(z % = B(z ! )Ui 

AOr 1 ) - 1 + Lri £(<r l ) = 2 hiz- { 

r=l 


where polynomials A and B are relatively prime, such that 
the horizon time for liftings have a total length N=m+l+p+ 
q+2n-l, and in which the number of elements o u in the 
liftings taken from the input signal is greater than the 
number of elements a y in the liftings taken from the output 
signal for output tracking in control systems and vice versa 
(cj u <g v ) for input tracking in communication systems. The 
extended horizon liftings for output tracking over a selected 
input signal window p u and a selected output signal window 
py are of the form 


r 



P«=[0 o, 1, 1, p. 





q 


n- 1 ^ 


1,1 1 , 0 0 , 0 , ..., 0 ] 


0,...,0, 0 0, 1,1 1J 


where m, 1, p, n, q, and n-1 are elements of the lifting, and 
m^0 and 1=^0 are arbitrary, q=0, p e R p is an arbitrary (or 
null) 0-1 vector chosen identically in both p M and p and n>0 
is the order of the irreducible plant. If the system is obtained 
by zero-order hold digitization of a continuous-time plant, 
the integer q^0 can be chosen arbitrarily. Since there are 
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6 


m-l-1 zeros on the left of p y and only m zeros on the left of 
p u , the number of elements c u is greater than a y for output 
tracking. For input tracking, the extended horizon liftings 
are of the same general form, but now there are m+1 zeros 
on the left of p M , and only m zeros on the left of p y so that 5 
the number of elements o u is less than a y . 

BRIEF DESCRIPTION OF THE DRAWINGS 


a logical “or” over the columns of S y , the mapping W is 
one-to-one. 

Using the above notation, the following “small” vectors 
are defined, 

Y,(k)^S y Y(k); S^W(p y )eR°> + " (8a) 

U,(k)^S K U(k); S„=W(pJeR°^ (8b) 


FIG. 1 graphically depicts partial horizon vectors U 5 and io 
Y s using windows p M and p r 
FIG. 2 is a block diagram of a generalized lifting system 
model G(p w ,p ;y ) defined by Equation 19 applicable to control 
systems in general. 

FIG. 3 illustrates graphically the determination of a small 15 
matrix H, from input and output signal liftings in respective 
windows p„ and p r 

FIG. 4 is a block diagram of a generalized lifting system 
model under zero annihilation conditions BS w r =0 and 
A(S/) r =0 defined by Eqs. (27) and (28). 20 

FIG. 5 is a block diagram of a zero annihilation periodic 
(ZAP) control law having a closed-loop system arising from 
the output tracking of extended horizon liftings of Equations 
(35) and (36). 25 

FIG. 6 is a pole-zero plot of a 12-state nonminimum phase 
plant with a sampling interval T=0.025sec in which poles are 
represented by “X” and zeros by “0.” 

FIG. 7 is a graph of open-loop output response of the plant 
in FIG. 6 to an initial condition. 30 

FIGS. 8a and 86 are graphs of respective output and input 
closed-loop response of the plant in FIG. 6 using deadbeat 
control with Lozano’s 2n-lifting (CT„=a y ) with a vibration 
damping objective. 

FIGS. 9a and 9b are graphs of respective output and input 35 
response using ZAP control of FIG. B with an extended 
horizon input lifting p u defined by m=0, 1=40, p=0, n=12, 
q=0, N=63. In this case a„=52 and cr y =12, so that G u >o y . 

FIG. 10 illustrates application of the present invention to 4Q 
communications for equalizing a nonminimum phase chan- 
nel. 


Y/(k)^S/Y(k); S/^W(p/)eR^ )+ * ( 8c) 

where pyand p u are specified 0-1 window vectors, and p y c 
is defined as the 0-1 complement of p r The vector Y/(k) in 
Eq. (8c) is denoted as the complementary output since it 
comprises all elements of the vector Y(k) which are not 
included in Y/k). 

A formula to reconstruct Y(k) from Y/k) and Y/(k) is 
now derived. First, it is to be noted that the quantity 
&=[S y r , (S/)T is a permutation matrix. Hence its inverse 
is its transpose, i.e., S J / 7 S ;y =I, which upon expanding gives 

S/S y HS y c ) T S y c =l (9) 

Multiplying both sides of Eq. (9) on the right by Y(k) gives 
the desired formula, 

Y{k)^ y T Y s {k)HS y c ) r Yf (*) (10) 

Example 2. FIG. 1 graphically depicts partial horizon vec- 
tors U s and Y s for the limited liftings case of n=3 and N=6, 
as determined by window vectors p u =[0, 1,1, 1,0,01 and p >; = 
[0,0,1,1,1,0]. In this case, p y c =[l,l, 0,0,0, 1] and one can 
compute, 


S U = W( p„) = 


0 1 0 0 0 0 

0 0 1 0 0 0 ; 

0 0 0 1 0 0 


Sy = W(p y ) = 


0 0 1 0 0 0 
0 0 0 1 0 0 
0 0 0 0 1 0 


DETAILED DESCRIPTION OF THE 

INVENTION 45 

A class of liftings (N>n) will first be defined by general- 
izing the lifting of Albertos (Eq. (4)). For this purpose, it will 
be useful to construct the “small” vector Y/k) from Y(k) as 
follows, 5Q 

Y^(k)=S > ,Y(k)eR <I >' 

where S y eR CTyX/v is a selection matrix which sifts out <j y 
elements of Y(k) for inclusion into Y^k). 

The matrix S y is most conveniently constmcted from a 55 
0-1 vector p y . For example, if p y =[l,l,0,0] then Y/k) e R 2 
contains the first two elements of Y(k) e R 4 and the selection 
matrix is given as, 


S y c = W( Py c ) = 


1 0 0 0 0 0 
oioooo 
oooooi 


In the present invention, a new family of extended liftings 
(N^2n-1) will be defined using the notation developed 
above for Albertos’ N-Hfting in Eq. (4). As a key step, it will 
be assumed that U(k) is chosen as zero outside the window 
defined by p u . Mathematically this can be written as, 

(I-S u T S u )U(k )= o (11) 

Consider the following nonminimal state- space realization 
of Albertos’ N-lifting Eq. (4) determined by using the small 
vectors in Eqs. (8) and identity Eq. (11), 


A systematic method to construct S y is defined as follows: 
form a diagonal matrix from the entries of vector p >( , and 
then remove all rows made up entirely of Os. This construe- 65 
tion defines the mapping WiR^^R^y*^ for which one can 
write S y =W(p y ). Since p y can be uniquely reconstructed by 


’ to i r a bsj i r u*“ i) i r m? 
urn J = [ 0 0 \ [u s (k~i) j + [ J 

r to i 


u s (k) 


( 12 ) 


(13) 


Substituting for the state in the output Eq. (13) gives the 
alternative output equation, 
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-continued 


[ Y(k-l) ] (14) 

Y s (k) = [SyA SyBS u T ) Us{lc _ l) +SyHS/U s (k) 

For analysis purposes, it is convenient to transform the 
open-loop plant using the similarity transformation T where, 



S/=W( p/) 


(15) 


10 

(16) 

(17) 15 


It can be verified that transformation T is square and 
invertible. Furthermore, T is in the form of a permutation 
matrix which reorders the state such that the components of 
appear first, and the remaining elements follow in the 20 
specified order. Since T is a permutation matrix, the inverse 
of T is given simply by its transpose (cf., Barnett Matrices: 
Methods and Applications, Clarendon Press, Oxford, 
England, 1990, pp. 374), i.e.. 


T' 1 =T 7 =[C t ,M t } 


(18) 


Transforming the open-loop dynamics of Eq. (12) by the 
similarity transformation of Eqs. (15) and (16), gives rise to 
a very useful representation denoted as the Generalized 30 
Lifting System Model, 


N 

r * "\ 

p y =[l,i,...,U] 

The 2n-lifting of Lozano-Leal, supra, is equivalent to the 
specific choice of periodic windows of horizon length N=2n 
as shown in Eq. (21), which is a special case also to be 
excluded from the present invention. 

N (21) 

n n ~ 1 

p M = [0,1,1 1, 0,...,0] 

n— 1 n 

Py — 10,0, . . . ,0, 1,1,..., 1] 

Since each choice of p and p v gives rise to a unique lifting, 
there are a total of 2 Ar+ " possible generalized liftings over a 
horizon of length N. Aside from the liftings of Eqs. (20) and 
(21), it appears that none of the other possible generalized 
liftings have been investigated in the literature. It will now 
be shown that many of the other possible generalized liftings 
have very useful properties, particularly those of N>2n-1. 

Zero Annihilation 

Focus will now be on those liftings for which the trans- 
mission zeros of the (squared down) lifted system lie at the 
origin. For notational convenience, we define the “small” 
matrix H v by, 

H-SyHSj ( 22 ) 


Generalized Lifting System Model, G(p u ,p y ): 


' Y s (k) 


SyASy T SyA (Sy C ) T 

SyBS/ 


' Ys(k- 1) 

Y s c (k) 

= 

Sy C ASy T S/A (S/) 7 ’ 

Sy C BS U T 


Y s c (k-l) 

U s (k) 


0 0 

0 


U 5 (k- 1) 


SyHSj 

S y c HS u r 

1 


40 

U s (k) 


It is noted that the generalized lifting G(p u , p y ) of horizon 45 
length N<n is defined uniquely by the choice of selection 
windows p„ and p >t , from which the matrices S u , S y ,S y c are 
calculated. The generalized lifting system model of Eq. (19) 
is depicted in the block diagram of FIG. 2. It is seen that 
and Y/ form two coupled subsystems which are driven by 50 
a common input Uy. It is also noted that the transmission 
zeros of the transfer function from U^(k) to Y^(k) are affected 
by the choice of windows p„ and p This is a key feature 
which will be used to advantage in applications of the 
present invention. 55 

The lifting G(p u , p y ) generalizes a number of existing 
results. For example, Albertos’ lifting in Eq. (4), taken from 
the paper of Albertos, supra, is equivalent to a choice of 
periodic windows (horizon length N) as shown in Eq. (20) 
for input and output liftings without regard to the plant order 60 
which should a priori be excluded for consideration in 
applications of the present invention. 


The quantity ¥L S appears in many expressions and will play 
an important role in subsequent proofs. The matrix can be 
obtained directly from p u and p y by writing p u along the top 
of H s and p y along the side of H. This arrangement is 
depicted graphically in FIG. 3. The matrix H s is then the 
submatrix defined by the elements of H having l’s along 
both borders. Two properties of interest concerning H, are, 
Output Tracking (OT) Condition: 




(23) 


Input Tracking (IT) Condition: 


H s r H=l (24) 

where I is the identity matrix and the superscript t denotes 
the Moore-Penrose inverse. It is noted that both the OT and 
IT conditions are satisfied if H^. is square and invertible. The 
OT condition is intended for gain adjustments in control 
systems, and the IT condition is intended for equalization of 
communication channels. 

The property of placing transmission zeros of the lifted 
plant to the origin is characterized in the following: 
Lemma 1 (Zero Annihilation) Assume that windows p u 
and p y satisfy, Zero Annihilation (ZA) Conditions: 

BS U T =Q (25) 

A(Sff = 0 (26) 

where, 


( 20 ) 


S U =W( pj; s~W( Py y, S/=W( p/) 


p u =tU U1 


65 Then, 

(i) the generalized lifting of Eq. (19) has a simplified 
representation, 
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Y s (k)^S y r Y s (k-l)+H s UM (27) 

Yf(k)=SfAS/Y s (k-l)+SfHS u r U s (k ) (28) 

as illustrated in FIG. 4, and 

(ii) Yf is unobservable from and has stable (deadbeat) 5 
dynamics. 

Furthermore, 

(iii) if Hy in Eq. (23) is square and invertible, then the 
transmission zeros of the lifted transfer function in Eq. (19) 
from Uy to are annihilated (i.e., lie at the origin), and 10 

(iv) if the OT condition Eq. (24) is satisfied then the 
transmission zeros of the “squared down” lifted transfer 
function Eq. (19) from V (where U^H/V) to Y^ are 
annihilated. 

Proof: Results (i) and (ii) follow by substituting the ZA 15 
conditions Eqs. (25) and (26) into Eq. (19) , to give, 


’ Y s (k) 


' Y s {k- 1) 

yjm 

=A P 

Y s c (k- 1) 

U s (k) 

L J 


U s (k-l) 


Y s (k- 1) 


Y s (k) = C p 


Ys c (k- 1) 

U s (k~ 1) 


+ D p U s (k) 


where. 



SyASy T 0 0 


SyHSj 

II 

Sy'AS/ 0 0 

£ 

II 

S/HS/ 


O 

O 

0 

1 


I 


C p = [S y AS y T 0 0]; D p = H s 


(29) 


20 

(30) 

25 

(31) 

30 


(32) 


According to standard definitions (cf., E. J. Davison and S. 
H. Wang, “Properties and calculation of transmission zeros 
of linear multivariable systems,” Automatica, Vol. 10, pp. 
643-658, 1974) values of X satisfying, 


det 


-U 



(33) 

40 


are the transmission zeros of the transfer function from Uy(k) 
to Y s (k) defined by the state-space model (A p ,B p ,C p t> p ). 
Consider the following identity, 


det(G)=det(G 22 )det(G n -G 12 G 22 1 G ! 2i) 


45 

(34) 


where, 




G11G12 

G21G22 


50 


and the inverse of G 22 is assumed to exist. Assuming that 
is invertible, the identity Eq. (34) can be applied to Eq. (33), 
which gives upon substituting Eqs. (31) and (32), 


det 


-X-I 

X 

X 


0 

-XI 

0 


0 

0 

-X- 1 


= (-X) NlG “ 


60 


which proves (iii). Result (iv) follows by an identical 
analysis assuming that the OT condition Eq. (23) holds, and 
that the simplified plant (A^B^C^D^,) in Eqs. (29)— (32) has 
been squared down by a precompensator H f. 

Lemma 1 is important since it gives conditions which p w 65 
and p y must satisfy for the generalized lifting G(p„, p y ) to 
have its transmission zeros at the origin. All results in 


10 

Lemma 1 can be simply understood by comparing the block 
diagram of FIG. 3 with FIG. 4 and noting all of the blocks 
that have vanished under the ZA conditions. It is seen that 
Yf no longer couples into the Y^ subsystem. Furthermore, 
the Yf subsystem has become deadbeat i.e., all of the poles 
of the Yf subsystem are at the origin. Most importantly, 
there is now only one forward path from Uy to Y^. Clearly 
something drastic has happened to the system zeros. Rigor- 
ously, if Hy is square invertible, result (iii) of Lemma 1 states 
that the transmission zeros of the transfer function from 
Uy(k) to Y 5 (k) have been placed to the origin (i.e., annihi- 
lated). If Hy is not square but the OT condition holds, result 
(iv) of Lemma 1 states that the zeros the lifted plant “squared 
down” by a precompensator H/ are annihilated. 

Extended Horizon Liftings 

The result of using extended horizons for the liftings 
introduces a new class of liftings which satisfy the condi- 
tions of Lemma 1. 

Theorem 1 A class of generalized liftings G(p u , p y ) which 
satisfy both the ZA and OT conditions is of the following 
form, 

Extended Horizon Lifting (OT Form): 




1 1. P> 




(35) 


(36) 


where elements mi?0 and 1^0 are arbitrary, q=0, p e R p is 
an arbitrary (or null) 0-1 vector chosen identically in both p u 
and p y ; and n>0 is the order of the irreducible plant of the 
form given by Eq. (1). Furthermore, if the plant (system) of 
Eq. (1) is obtained by a zero-order hold (ZOH) digitization 
of a continuous -time system, the integer q^0 can be chosen 
arbitrarily. 

Proof: By definition, the matrix B in Eq. (5c ) has the 
sparse form B=[O fc !X b ] where O z ,eR ArxAr-n+1 is a matrix of all 
“0” elements, and X^eR^*" -1 . By construction of S u from p„ 
in Eq. (35), the nonzero elements of S„ multiply only 
elements of 0^ in the product BS u r . Hence BS/^0. Like- 
wise, the matrix A in Eq (5a ) has the sparse form A=[OJX J 
where O^eR^^ - ” is a matrix of all “0” elements, and 
X a eR jVxn . By construction of S y c from p y c (i.e., 0-1 comple- 
ment of p y in Eq. (36)), the nonzero elements of S v c multiply 
only elements of O a in the product A(S >( c ) r - Hence A(S 3 , c ) r =0 
and the ZA conditions of Eqs. (25) and (26) are satisfied. 

In order to show the OT condition, first consider the case 
where q^0. Then choice of Eqs. (35) and (36) ensures that 
the small matrix Hy is of the form shown in FIG. 3, 



where 3 e R p>9? is lower triangular with a nonzero diagonal 
(and hence is invertible), and Ke R" xn is given by. 
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( 38 ) 


<H = 


5 


^ hln+q-\ ■ • • hn+q J 

Since 3 is full rank, it follows from the special structure of 
Eq. (37) that is full rank if has full rank. To show that 
tK has full rank, let (A,b,c) be any minimal (i.e., controllable 10 
and observable) state-space realization of the transfer func- 
tion Eq. (1). The Markov parameters {h,} can be written as 
h=cA l_1 b,i=l, . . . , «. Substituting into Eq. (38) gives. 


cA n ^b 


H = 


cA q b 


[ cA^^b . . . cA n+q ~ ] b 
M = OA q di 


( 39 ) 


15 


20 

( 40 ) 


where I is a reversed identity (i.e., I=[e l9 . . . ,e„], I=[e n , . . 

. ,ej ]), and O and C are observability and controllability 
matrices of (A,b,c), respectively. Since (A,b,c) is control- 25 
lable and observable, it follows that O and C are each full 
rank. The fact the system or plant of Eq. (1) is obtained by 
a ZOH digitization implies that AeR nx " is full rank (i.e., A 
is a state-transition matrix). These facts together imply that 

30 

‘K in Eq. (37) is full rank and hence is invertible. 

If the system of Eq. (1) is not obtained by ZOH digiti- 
zation, the matrix A may not be full rank. However, the 
results still hold with the restriction that q=0 since in this 
case A q is replaced by identity matrix I in relation Eq. (40), 
which is always full rank. 35 

Theorem 1 is important because it extends Lozano’s 
2n-lifting to a much larger class of liftings which enjoy the 
same zero annihilation properties, and thus provides greater 
latitude (extra degrees of freedom) in designing control 
systems. Note that the extended horizon liftings have total 40 
horizon length N=m+l+p+q+2n-l, which can be readily 
chosen longer than Lozano’s lifting for which N=2n. Hence 
the phrase “extended liftings horizon”. It will be seen that 
the extra degrees of freedom provided by N=k+2n-l , where 
k=m+l+p+q, for example, will overcome several difficulties 45 
associated with the 2n-lifting. For convenience, this new 
class of liftings is referred to herein as having N^2n-1. 

The following result is essentially a “dual” to the previous 
theorem. 

Theorem 2 A class of generalized liftings G(p tt , p v ) which 50 
satisfy both the ZA and IT conditions is of the following 
form, 


Extended Horizon Lifting (IT Form): 


N 




m Ip n q n- 1 

p u =[0,...,0, 0,...,0, p, 1,1 1, 0,...,0, 0,...,0j 


55 


( 41 ) 


60 


l p q n - 1 

p y =[0,...,0, ~p? oT^7,0, 0 0, 1,1, 1] 


( 42 ) 


where rn^O and 1^0 are arbitrary, q=0, peR^ is an arbitrary 65 
(or null) 0-1 vector chosen identically in both p„ and p v ; and 
n>0 is the order of the irreducible plant Eq. (1). Furthermore, 
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if the plant or system Eq. (1) is obtained by a zero-order hold 
(ZOH) digitization of a continuous-time system, the integer 
q^0 can be chosen arbitrarily. 

Preof: The ZA conditions follow using the same argu- 
ments found in the proof of Theorem 1. Consider the IT 
condition Eq. (24). Given the lifting Eq. (41), H, has the 
form, 


p n 


( 43 ) 


p 

H s = n 
/ 


3 

X 

X 



From the structure of Eq. (43) it follows that H^ has full rank 
if both 3 and Khave full rank. The remainder of the proof 
is identical to the proof of Theorem 1. 

Note: Properties of Lozano’s 2n-lifting (N=2n) follows 
directly from its interpretation as a special case of the liftings 
N^2n-1 in Theorem 1 and Theorem 2. 

Corollary 1 Lozano’s 2n-lifting Eq.(21) satisfies the ZA, 
OT and IT conditions. 

Proof: The lifting Eq. (21) is equivalent to the special case 
of the liftings in Theorems 1 and 2 where m=0, 1=0, p=l , 
P=[0], q=0. 


ADVANTAGES 

The advantages of the extended horizon liftings in Eqs. 
(35) and (36), and in Eqs. (41) and (42) over Lozano’s 
2n-lifting are as follows: 

(i) If one chooses 1>0 in Eqs. (35) and (36), there are more 
control inputs than outputs in the lifted system (i.e., G„>G y ). 
It is shown below with respect to application to plant inverse 
control that these extra degrees of freedom can be used to 
design a controller which minimizes a quadratic control cost 
while simultaneously satisfying a deadbeat tracking objec- 
tive. This significantly reduces control gains compared with 
Lozano’s lifting. 

(ii) If one chooses 1>0 in Eqs. (41) and (42), there are 
more outputs than inputs in the lifted system (i.e., 0y>0„). It 
is shown below with respect to channel equalization in 
applications to communications that these extra degrees of 
freedom can be used to minimize a quadratic error when 
estimating the input from measurements of the output. This 
is significant for reducing noise in problems of nonminimum 
phase channel equalization. Furthermore, the p vector can 
contain additional message information to increase the chan- 
nel throughput. 

(iii) If one chooses m>0 in Eqs. (35) and (36) or Eq. (37), 
there is an extra m*T seconds of free time which can be used 
to perform computations (where T is the sampling interval). 
Since m can be chosen arbitrarily, the use of extended 
horizon liftings for either control or equalization applica- 
tions is not constrained by real-time computer limitations. 
This is particularly useful for adaptive implementations 
which involve additional computation. 

APPLICATION TO PLANT INVERSE CONTROL 

The placement of the transmission zeros to the origin by 
the class of extended horizon liftings of the form given by 
Eqs. (35) and (36) allows stable invertibility of the transfer 
function from U J (k) to Y 5 (k). A control law which will be 
discussed next deadbeats the response Y s (k) to follow the 
desired Y d (k), subject to the minimization of a quadratic 
control cost. 

To derive the desired controller, define the output error as, 



E(k)=Y d (k)-Y s (k) 
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-continued 


Substituting Eq. (27) into Eq. (44) gives, 
E(k)=~SyAS y r Y s (k~l)-H s U s (k)+Y d (k) 


(44) 


(45) 


Bel - 


H S H S ‘ 

SfHS/HJ 

HJ 


Consider the problem of forcing the error in Eq. (38) to zero 
in a single step, while minimizing a quadratic control cost 
penalty, i.e., 

(46) 10 

min U/(k)U s (k) 

U s {k ) 

subject to 

E(fc)=0. (47) 15 


In light of the output tracking condition of Eq. (23), this 
minimization problem can be solved uniquely (see D. S. 
Bayard and D. Boussalis, “Noncolocated structural vibration 
suppression using Zero Annihilation Periodic control,” 2nd 
IEEE Conference on Control Applications, Vancouver 
Canada, September 13-16, 1993) to give: 

Zero Annihilation Periodic (ZAP) Control Law: 


Uf^Hji-SyAS/Y^k-iy+Y^k)) 

=K°Y s (k-l)+L°YM 


where the corresponding feedback gains are defined as, 

^-H/SyASy 7 ' 

L°=HJ 


‘ Y s (k) 


’ Y s {k- 1 ) 

Y s c (k) 

—Ad 

Y s c (k- 1 ) 

_ U s (k) _ 


_ UAk-D 


where, 


A c i - 


Sy C (I ~ HS U T H s , Sy)ASy 1 
-H/SyASy 7 


+ BdY d {k) 


0 0 
0 0 
0 0 


20 


(48) 25 

(49) 


(50) 

(51) 
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Here the superscript “o” is chosen to emphasize the fact that 
the control nulls (i.e., deadbeats) the output. Also, in light of 
the OT condition, has full row rank and one can write 
Z)- 1 (cf., Barnett, supra ). 

For convenience the ZAP control law is summarized in 
the block diagram of FIG. 5. It provides the following result. 

Lemma 2 (ZAP Control) Consider the closed-loop system 
arising from the OT extended horizon lifting Eqs. (35) and 
(36) under ZAP control Eq. (49). Then, 

(i) The quadratic control cost Eq. (46) is minimized at 
each stage, subject to the deadbeat tracking constraint Eq. 
(47), 

(ii) All closed-loop poles are at the origin (i.e., the 
closed-loop response is deadbeat), and hence Y s (k) con- 
verges to Y d (k) in a single step. 

(iii) The closed-loop system is internally stable (e.g., 
Y/(k) remains bounded). 

Proof: Result (i) follows from the OT condition in Eq. 
(23) and well known minimum-norm properties of the 
Moore-Penrose inverse (cf., Barnett, supra). Now form the 
closed-loop system from the simplified lifted plant Eqs. (27) 
and (28) under ZAP control Eq. (49), 


35 


40 


45 


50 


55 


(52) 


60 


(53) 


65 


Results (ii) and (iii) follow by substituting the OT con- 
dition Eq. (23) into Eq. (53), and noting that the resulting 
closed-loop matrix is stable with all of its eigen values at the 
origin. 

Result (i) of Lemma 2 is important because it indicates 
that control gains associated with using extended horizon 
liftings of the class denoted herein as N^2n-1 will be 
significantly reduced compared to those from using Loza- 
no’s lifting (N=2n). Result (iii) of Lemma 2 is important 
because it ensures that the complementary output Y/ 
remains “well behaved” even though it is not being con- 
trolled directly. 

Remark 1 Instead of deadbeat control, a pole placement 
scheme can be obtained by modifying the deadbeat con- 
straint Eq. (47) to become E(k)=aE(k-l) in which case the 
ZAP control becomes U/(k)=K°Y,(k-l)+L°Y rf (k)-~ 
otL°E(k-l) 

Example 3 As an example, a 12-state nonminimum phase 
transfer function is shown in FIG. 6 (pole-zero plot, sam- 
pling time T=0. 025 sec), adopted from the ASTREX flexible 
structure model(cf. Bayard, 1993, supra and A. Das, J. L. 
Berg, G. A. Norris, D. F. Cossey, T. J. Strange IB, and W. T. 
Schlaedgel, “ASTREX — A Unique test bed for CSI 
research,” Proc. 29th IEEE Conference on Decision and 
Control, Honolulu, Hawaii, pp. 2018-2023, December 
1990). The open-loop response to an initial condition is 
shown in FIG. 7. A simulation is first run using the deadbeat 
control Eq. (49) with the Lozano’s 2n-lifting, m=0, 1=0, p=l, 
p=[0], n=12, q=0, N=2n=24, and using Y^0 (i.e., a vibra- 
tion damping objective). The response is shown in FIG. 8 to 
reach -3x1 0 4 at the output and 500 at the input. As expected 
from the theory, the response is deadbeat after a single 
horizon. However, this control law is unusable since the 
units are in volts, and the allowable range is only ±10 Volts. 
It is emphasized that these extraordinarily large responses 
are typical of the 2n-lifting due to the fact that the plant is 
inverted on a horizon of length N.T=0.6 sec. 

In an attempt to get a practical response with the 2n-lift- 
ing, the pole placement control of Remark 1 is used. The 
choice a=0.5 is made to get approximately a 10-second 
decay time (any slower would be worse than the open-loop 
response). The results are simulated but not shown here 
since it turns out that the responses are reduced 50%, and are 
still unacceptably large by several orders of magnitude. The 
ZAP control using an extended horizon lifting of the class 
N^2n-1 with m=0, 1=40, p=0, q=0, n=12, N=63 was tried 
next. The results are shown in FIG. 9 where it is seen that 
both the input and output are well within the allowable 
ranges. The deadbeat nature of the response is also noted, as 
the vibrations are damped instantaneously after the first 
horizon at time T=3sec. This demonstrates the advantage of 
using the extended horizon liftings of the present invention 
instead of Lozano’s 2n-liftings. 

Application to Channel Equalization 

An important problem in communications is that of 
equalizing (i.e., compensating for undesired amplitude-fre- 
quency response) a nonminimum phase channel. The usual 
problem is that the channel cannot be inverted in a stable 
fashion. However, using the input tracking (IT) extended 
horizon lifting form of Eqs. (41) and (42), this problem can 
be overcome. 

Referring to FIG. 10 , let U 5 (k) be the sequence of mes- 
sages to be sent by a transmitter T, and assume that an IT 
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extended horizon lifting form given by Eqs. (41) and (42) is 
used to transmit the data (i.e., the signal sent is given by 
U(k)=S u r U :s .(k)). The channel C is assumed to be a stable 
linear nonminimum phase transfer function for a plant of the 
form given by Eq. (1) with a plant order n. Then from 5 
Theorem 2 and Lemma 1, the plant (channel) dynamics 
given by Eqs. (27) and (28) are 


Y x (k)=A s Y s (k-\)+H s U s (k) 

where 



A=SyAS y T * (55) 

At the receiver R, the quantity Y^(k)=S y Y(k) is measured, 
and it is desired to estimate the messages U^(k) which were 15 
sent. For this purpose, an output prediction is Y s (k) is 
formed as, 


? s (k)=A s Y s (k-l)+H s U s tt) (56) 

20 

and an estimate \J S found by minimizing the least squares 
criteria. 


(57) 

rain (Kj- Ys) t (Y s -Y s ) 25 

U s 

Since the IT condition Eq. (24) holds for the extended 
horizon lifting, the unique solution to Eq. (56) is given by, 

^^(A^.CA-l)-^^)) (58) 30 

where one can write H/=(H/H 5 ) -1 H^ r since Hy has full 
column rank (cf., Barnett, supra). It is noted that this 
dynamical system represents a stable inversion of the non- 
minimum phase channel characteristics. Furthermore, the 35 
use of an extended horizon lifting with 1>0 has allowed 
channel inversion by least squares which provides smooth- 
ing in case of noise. Lozano’s 2n-lifting in this application 
would not allow smoothing and would be very susceptible to 
noise. The use of p in the IT extended horizon lifting also 4Q 
provides a means for sending additional information with 
each packet, to improve the overall throughput. 
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noise. It is worth noting that this channel equalization 
approach can be made adaptive by using standard recursive 
algorithms since the new liftings are linear-in-the-param- 
eters. 

It is expected that the present invention (extended horizon 
liftings) will be useful in many areas of modem control, 
neural control, fuzzy control, adaptive control, communica- 
tions, adaptive filtering, signal processing, or other applica- 
tions where a stable system inverse is desired but not 
possible due to nonminimum phase constraints. 

I claim: 

1. An improved method for providing periodic gain 
adjustment in a plant of irreducible order, such as for control 
of electrical, mechanical or chemical systems, or for equal- 
ization of communication channels, comprises the choosing 
of horizon time N>n of liftings in periodic input and output 
signal windows p„ and p y , respectively, of said plant where 
N is an integer chosen to define the extent of each of said 
windows p u and p y and n is the order of the irreducible 
input/output plant model given by an equation of the form 

A(z-% ( = flOf')K ( 

A(z- J ) = l+ X at-*-, B(z~ l ) = ibtz- 1 
i=l r=l 

where polynomials A and B are relatively prime, such that 
the horizon time for liftings have a total length N=m+l+p+ 
q+2n-l, and in which the number of elements o u in the 
liftings taken from the input signal is greater than the 
number of elements g v in the liftings taken from the output 
signal for output tracking in control systems, and in which 
the number a u in the liftings taken from the input signal is 
less than the number of elements <j y taken from the output 
signals for input tracking in communication systems, and the 
following conditions hold for output tracking: 
where the superscript t denotes the Moor e-Penrose inverse 
of a submatrix U s obtained directly from p u and p y by 
placing said vector p H along the top of a matrix H and p y 
along the side of said matrix H and taking as said submatrix 
Hy all of the elements of said matrix H along both borders, 
and 

BSJ = o 


CONCLUSIONS A(S y c ) r =Q 


A general class of liftings N>n has been denominated as 
Ni?2n-1 in order to distinguish from the special case of 
Lozano (N=2n) and have been shown to have the same 
desirable zero annihilation properties of Lozano’s 2n-lifting. 

In contrast to Lozano’s 2n-lifting, all the new liftings include 50 
horizons greater than 2n, i.e., they are of the extended 
horizon type N>2n-1 and g u >g v or g^g,,. The use of 
extended horizons of the type N>2n-1 resolves many dif- 
ficulties associated with Lozano’s 2n-lifting. For example, a 
Zero Annihilation Periodic (ZAP) controller is defined for 55 
which the control gains can be significantly reduced as 
compared to Lozano’s lifting. This is due to a quadratic 
control cost which is minimized simultaneously with the 
deadbeat tracking objective. The effectiveness was shown in 
a simulation example where the control torque was reduced 60 
4 orders of magnitude. 

As a dual result, it was shown that a related class of 
liftings enables the on-line equalization of nonminimum 
phase channels in communication systems. This overcomes 
the standard bottleneck of inverting the channel in a stable 65 
fashion. Here, channel inversion is accomplished by least 
squares estimation which provides smoothing in the case of 


where S„-W(p w ); S y =W(p y ); S/=W(p v c ) and said liftings 
have the following representation 

y 5 (/t)=S 3 AS/k(*- 1 )+H s U s {k), and 

Y/(k)=S y c AS y T Y s {k-\)+S y c HS/U s (k ) 

whereby the transmission of zeros of said liftings from U v to 
Y s are annihilated and the transmission of zeros from V, 
where H/=H v V, to Y s are annihilated and H y =S >r HS„ r , and 
the following conditions hold for input tracking: 

H S T H~I 

BSj = 0 

A ( S /) 7 -0 

and said liftings have the following representation: 
Y s (k)^yAS y r Y s (k-\)±S/HSjU s (k\ 

and controlling the performance of said plant in response to 
said chosen horizon time. 
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2. An improved method as defined in claim 1 wherein 
extended horizon liftings for output tracking over a selected 
input signal window p u and a selected output signal window 
p y are of the form 


N 




m l p n 

p« = [0, . . . ,0, P, 1,1,..., 1, 



n - 1 



5 


10 


m 


Py- [0, . . . ,0, 


l p q 

0, . . . ,0, p, 0, . . . ,0, 




U/(k)=K°Y s (k-l)+L 0 Y d (k)-aL 0 E(k-\ ) 


where, 


E(k-D=Y d {k)-Ys(k) 

K°=-H s 'S y ASy T 

L°=HJ 

and a is a scalar less than unity. 

5. An improved method as defined in claim 1 wherein said 
extended horizon liftings are of the form, 


where m, 1, p, n, q, and n-1 are elements of the lifting, and 15 
m^0 and 1^0 are arbitrary, peR p is an arbitrary or null 0-1 
vector chosen identically in both p„ and p v and n>0 is the 
order of the irreducible plant, and q=0, unless the plant 
liftings are obtained by a zero-order-hold digitization of a 
continuous-time analog signal, in which case q>0 can be 20 
chosen arbitrarily. 

3. An improved method as defined in claim 2 wherein 
output tracking for control is implemented using the follow- 
ing Zero annihilation Periodic control law 
6 25 

Uf(k)=Hj(-S^S y r Y s (k-l)+Y^k)) 
where the feedback gains K° and If are defined as, 

K°--HjSyAS y r 30 




N 



-N 


n-1 


0,...,0, p, 1,1,..., 1, 0,...,0, 0,...,0] 

l p q n-1 n 

£7^1, "-p? oT^To, oTTTa 1,1, 7 . j] 


for input tracking such that there are m+1 zeros on the left in 
the window p u , and only m zeros on the left in the window 
p v so that the number of elements in the window o u is less 
than in the window a y . 

6. An improved method as defined in claim 5 for input 
tracking wherein an estimate of the transmitted word \J S is 
computed at the receiver of a communication channel as 


U>H 1 

U s =-HfA s Y s (k-l)-Y s (k)) 

where H/ is the Moore-Penrose inverse of H^, and Y d is a 

desired output signal to be tracked. where H/ denotes the Moore-Penrose inverse of H^. 

4. An improved method as defined in claim 2 wherein 35 
control is implemented using a pole placement scheme, ***** 



